Synthesis and Characterization of Single Crystal Samples of Spin-| Kagome Lattice 
Antiferromagnets in the Zn-Paratacamite Family Zn a .Cu 4 _ x .(OH) 6 Cl 2 
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The Zn-paratacamite family, ZnaCxu-x (OHQsCh for x > 0.33, is an ideal system for studying 
spin-| frustrated magnetism in the form of antiferromagnetic Cu 2+ kagome planes. Here we report 
a new synthesis method by which high quality millimeter-sized single crystals of Zn-paratacamite 
have been produced. These crystals have been characterized by metal analysis, x-ray diffraction, 
neutron diffraction, and thermodynamic measurements. The x = 1 member of the series displays 
a magnetic susceptibility that is slightly anisotropic at high temperatures with \c > Xat- Neutron 
and synchrotron x-ray diffraction experiments confirm the quality of these x — 1 single crystals and 
indicate no obvious structural transition down to temperatures of T = 2 K. 

PACS numbers: 81.10.-h, 75.50.Ee, 61.05.cp 



Geometrically frustrated magnetism 1 ^ is a forefront 
area of research in condensed matter physics, as such sys- 
tems offer a unique terrain in which to search for novel 
magnetic ground states. The spin-g nearest- neighbor 
Heisenberg antiferromagnet on the kagome lattice, which 
consists of corner sharing triangles, is a particulary 
promising system in which to search for unique quantum 
phases including the "resonating valence bond" (RVB) 
state proposed by Anderson 3 - or other quantum spin liq- 
uid states. A broad theoretical and numerical consen- 
sus has emerged that the ground state of this system is 
not magnetically ordered 4 - - — , with a variety of proposed 
ground states including gapped spin liquids^, gapless 
spin liquids*^, and valence bond solid (VBS) state o 11 ' 12 . 
However, experimental investigation on this system has 
long been hampered by the fact that most early realiza- 
tions of the kagome lattice antiferromagnet feature either 
large spins or structural distortions. 

The material ZnCu 3 (OH)gCl 2 13 ' 14 is among the best 
realizations of a spin-i kagome lattice antiferromagnet 
yet synthesized. This material is a member of the Zn- 
paratacamite family Zn^Cu^^OH^C^ with x = 1. 
With lattice parameters a = b = 6.83 A and c = 14.05 A, 
ZnCu3(OH)gCl2 is rhombohedral (trigonal setting) and 
consists of kagome lattice planes of spin- h Cu 2+ ions sep- 
arated by layers of non-magnetic Zn 2+ ions as shown in 
Fig. 1(a). Measurements on powder samples 1 ^ have found 
no sign of long range order or spin freezing down to tem- 
peratures of 50 mKi£r— , despite a strong antiferromag- 
netic superexchange interaction of J w 17 meV and a 
Curie- Weiss temperature of &cw = —300 ± 20 K. There 
is no evidence of a spin gap down to at least J/200 16 ' 19 ' 20 . 
The most significant deviations of this material from 
the idealized model are likely the presence of about 5% 
weakly-coupled Cu 2+ ions lying on out-of-planc metallic 
sites, which may be responsible for the Curie-like suscep- 
tibility at low temperatures, and Dzyaloshinksii-Moriya 
or exchange anisotropy interactions. Recent anomalous 
x-ray diffraction measurements indicate that dilution of 




FIG. 1: (color online) (a) Structure of ZnCu 3 (OH) 6 Cl 2 with 
only Cu 2+ (large brown spheres) and Zn 2+ (small red spheres) 
displayed. The Cu-Cu bonds (thick black solid lines) are 
all equivalent as are the Cu-Zn bonds (thin green dotted 
lines), (b) A single crystal sample of ZnCu3(OH)6Cl2- (c) 
A schematic of three zone furnaces. The red bars indicate the 
positions of the thermocouples. 



the kagome plane sites with Zn ions is not significant^!. 
Studies on Zn-paratacamite sample o 17 ' 22 ' 23 with x < 1 
can be useful in determining the effect of any out-of-planc 
impurities that might be present even in the nominal 
x = 1 samples, as well as serving as model systems for 
weakly coupled kagome planes. 

Previous studies of Zn-paratacamite have utilized pow- 
der samples grown by hydrothcrmal methods. The fail- 
ure of these methods to produce large single crystals has 
been ascribed partly to the low decomposition tempera- 
ture of ZnCu3(OH)eCl2. However, the primary reason for 
the lack of single crystal synthesis is likely the fact that 
the reported synthesis 1 ^ produces bubbles of CO2 , result- 
ing in an unstable crystallization environment. Further 
understanding of the spin behavior of ZnCu3(OH)eCl2, 
and further insight into the ground state of the spin-i 
kagome lattice antiferromagnet, will require studies on 
single crystal samples. Here we report a new synthe- 
sis method by which high quality millimeter-sized single 
crystals have been successfully produced. These samples 
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have been characterized by a variety of measurements. 

Single crystal samples of Zn-paratacamite, 
Zn x Cu4_ x (OH)gCl2, were grown hydrothcrmally in 
furnaces that were setup similarly to those used to grow 
small single crystal samples of the atacamite family^ 
and Mg 2; Cu4_ 2 ;(0H)6Cl2^. Here, starting materials of 
CuO, ZnCl2, and H2O, in amounts listed in Table HI 
were charged into a fused quartz tube (ID 6 mm, OD 
13 mm for x=0.8 and 1.0 or ID 9 mm, OD 15 mm 
for x=0.9). The quartz tube was sealed after purging 
air with a mechanical pump. The sealed quartz tube 
was prereacted for two days in a box furnace at 185 
°C. After prereaction, a green-blue microcrystallinc 
powder was formed. Powder x-ray diffraction measure- 
ments of this product have indicated the presence of 
Zn^Cu^^OH^C^. This shows successful synthesis by 
the reaction 

(4-.t)CuO + ZnCl 2 + 3H 2 -> 

Z n:r Cu4_ x (OH) 6 Cl2 + (1 - x)ZnO. (1) 

This synthesis without the production of CO2 suggests 
the possibility of an environment stable enough for single 
crystal growth. 

Millimeter-sized single crystals were synthesized 
through a recrystallization process in a three-zone gradi- 
ent tube furnace. A schematic of such a furnace is shown 
in Fig. (He). The sealed, prereacted quartz tubes were 
placed horizontally into the furnace at room temperature. 
The furnace temperature was isotropically increased to a 
fixed temperature, ranging from 165 °C to 180 °C in var- 
ious reactions. The temperature of the cold end was then 
slowly lowered. The sample and all growth parameters 
were undisturbed for roughly 20 weeks until large crys- 
tals were formed at the cold end. In the region where 
the crystals nucleated and grew, the temperature gradi- 
ent was measured to be approximately 1 °C/cm. At the 
end of the synthesis, the sample tubes were cooled down 
to room temperature at 1 °C/min. Crystals were then 
rinsed with deionized water, dried in air and kept in a des- 
iccator for storage. No decomposition of the crystals has 
been observed in air, water, or acetone. Precise control 
of the starting concentrations of CuO and ZnCl2 allows 
for synthesis of samples with variable Zn concentration, 
x. Data on crystals with x = 0.8, 0.9, and 1.0 are shown 
in Table HI More syntheses than listed in Table U were 
performed with ZnCl2 to CuO molar ratios ranging from 
2 to 10 and with ZnCl2 to H2O concentrations ranging 
from 1.2 mmol/ml to 7.7 mmol/ml. However, the x val- 
ues of the final products were fairly stable over this range 
of starting concentrations. At a fixed ZnCl2 to CuO ra- 
tio, the x value of the product increased with increasing 
ZnCl2 concentration, from an x value of 0.8 with a ZnCl2 
concentration of 1.2 mmol/ml to an x value of 1.0 with a 
ZnCl2 concentration of 2.8 mmol/ml. Products with an 
x value of 1.0 were also obtained for starting Z11CI2 con- 
centrations up to 5.6 mmol/ml, while even higher Z11CI2 
concentrations resulted in a slight lowering of x. This 
diminishment of x with very high ZnCl2 concentrations 
is likely due to a more acidic pH in those reactions which 
dissolves more CuO. The ZnCl2 to CuO molar ratio of 
the starting products had no obvious effect on the x val- 
ues of the product over the range of syntheses performed. 
The compositions of the crystals were measured by metal 



TABLE I: Growth and crystallography data. All samples have 
rhombohedral crystal system in R3m space group (a=/3— 90° , 
7=120°. Single crystal x-ray diffraction was performed at 
T — 100 K and refined by full-matrix least-squares on F 2 
with goodness-of-fit (GOF) listed. Curie- Weiss temperatures 
were determined from high temperature susceptibilities. 





x = 0.8 


x = 0.9 


x = 1.0 


Starting 


CuO (0.130) 


CuO (0.346) 


CuO (0.235) 


materials (g) 


ZnCl 2 (0.686) 


ZnCl 2 (2.985) 


Z11CI2 (2.015) 




H2O (4.0 ml) 


H2O (10.0 ml) 


H2O (4.5 ml) 


Hot zone temp. 


165 °C 


165 °C 


180 °C 


a 


6.8300(13) A 


6.8345(9) A 


6.8332(12) A 


b 


6.8300(13) A 


6.8345(9) A 


6.8332(12) A 


c 


14.029(3) A 


14.0538(19) A 


14.066(2) A 


Volume (A 3 ) 


566.77(19) 


568.51(13) 


568.80(17) 


Density, p 


3.765 g/cm 3 


3.759 g/cm 3 


3.755 g/cm 3 


total reflections 


2673 


3718 


5908 


indep reflections 


238 


225 


240 


GOF 


1.275 


1.264 


1.237 


Qcw 


-266(10) K 


-290(10) K 


-296(10) K 



analysis taken with an inductively coupled plasma atomic 
emission spectrometer (ICP-AES) with an error of ±0.04 
on x. Standards were prepared from commercially pur- 
chased solutions from Sigma-Aldrich, specific for ICP- 
AES measurements and designated as Trace SELECT 
grade or better. Five to ten well rinsed small single crys- 
tals from each synthesis tube, approximately 0.1 mg each, 
were dissolved into 2% w/w nitric acid for measurement. 
Unlike powder samples, the ease to rinse single crystals 
dramatically reduced the ambiguity from possible chem- 
ical contamination. The relative amounts of Cu and Zn 
determined from ICP metal analysis were used to calcu- 
late the values of x listed in Table |TJ Our nominal x = 1 
sample has previously been determined to have structural 
composition (Zno.85Cuo.i5)Cu3(OH)eCl2 via anomalous 
x-ray scattering^ Here, the kagome planes are fully oc- 
cupied with Cu, and anti-site disorde r 26 ' 27 with Zn on 
the Cu kagome site is not apparent. 

Single crystal x-ray diffraction was performed on a 
three-circle diffractometer coupled to a CCD detector. 
All samples were refined in the rhombohedral space group 
R3m (trigonal setting) and with lattice constants consis- 
tent with previous reports. The largest (1 to 10 mm 3 ) 
crystals typically grow as a bar-shape similar to the crys- 
tal shown in Fig. 1(b). From x-ray diffraction on more 
than 30 crystals, all of the side long faces were indexed 
as (1 1) (the faces were normal to the (1 1) recip- 
rocal lattice vector). Some of the smaller (< 0.5 mm 3 ) 
crystals were octahedrally shaped, with all eight faces 
indexed as (1 1). Based on these observations, we 
propose the following growth process: during the early 
stage of crystallization, primitive micrometer-sized crys- 
tals form as twelve-faced polyhedra due to the symmetry 
of the (1 1) directions. As the crystals grow larger in 
size, eight of the twelve (1 1) faces, possibly due to the 
specific local hydrodynamic condition, grow faster which 
transforms the crystals into larger, sub-millimeter octa- 
hedra. Eventually, four of these eight (10 1) faces form 
the four large side faces of millimeter-sized bar-shaped 
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crystals. 

The single crystal samples obtained by this synthe- 
sis were characterized by a variety of methods. Single 
crystal susceptibility measurements of an x — 1 sample 
were performed on a SQUID magnetometer (Quantum 
Design) using a 55.5 mg sample (different from the one 
in Fig. [T](b)) with an almost cubic shape (2.3 mm x 2.5 
mm x 2.7 mm). Figure. [U(a) shows the bulk suscepti- 
bility, x = M/H, as a function of temperature. The 
susceptibility of the crystal closely follows that measured 
on a powder sample with x — 1 (also plotted), which 
confirms that the stoichiometry and homogeneity of the 
single crystals match that of previously measured pow- 
ders. The inset of panel (a) shows the specific heat of an 
x = 1 single crystal measured at (j,qH = T and com- 
pared with data from a powder sample. Again, the two 
curves show good agreement, where the higher quality 
of single crystal measurements is likely due to the better 
thermal contact with the sample holder. In Fig.^b), we 
show the magnetic susceptibility plotted as MT/H, with 
an applied magnetic field of (J-qH = 1 T aligned along the 
different crystallographic axes, (x c ) and(xa&)- Powder 
data are also plotted (xpowder) for comparison. The cal- 
culated powder averaged value of the single crystal data 
(Xav = §Xa6 + \Xc) is indicated by the solid line. As ex- 
pected, the calculated Xav and Xpowder show good agree- 
ment. At high temperatures, the single crystal data show 
clear magnetic anisotropy with Xc > Xab- The inset to 
panel (b) shows the high temperature (80 K and 300 K) 
magnetization; the anisotropy is field independent, with 
M c /M a t, a constant up to /in-ff = 5 T. This high tem- 
perature anisotropy is consistent with the qualitative re- 
sults obtained from powder samples that were partially 
oriented in a magnetic field^ and /iSR measurement on 
single crystal samples^. The anisotropy of the magneti- 
zation should play a role in elucidating the effects that the 
Dzyaloshinskii-Moriya interaction, exchange anisotropy, 
or out-of-plane impurities have on the low temperature 
physics of the materia l 30 ' 31 . The bulk susceptibility of 
crystals with x = 0.8 and 0.9 are also consistent with 
powder samples (not shown). A detailed study of the 
susceptibility and specific heat of samples with different 
x- values is a topic of ongoing investigation. 

Among the various predicted ground states for the 
spin-i ncarest-neightbor Heisenberg kagome antiferro- 
magnet, several theoretical studies have suggested a 36- 
site valence bond soli d 11 ' 12 . This proposed ground state 
features a s/12 x s/12 enlargement of the unit cell, with 
the 36 spin sites paired into 18 nearest-neighbor dimcr 
singlets, where 6 dimers lie around a central pinwheel 
configuration while another 6 lie around hexagons. The 
two possible coverings of the pinwheel are degenerate to 
high order, while sets of three dimers can resonate around 
a perfect hexagon^. The other 6 dimers are presumed 
to be static. It has been suggested 3 ^ that a VBS order 
might lead to a slight structural distortion in which the 
distance between two magnetic ions paired in a singlet is 
diminished. Such a distortion, if large enough, could in 
principle be measured in a synchrotron x-ray experiment. 

We performed diffraction measurements on our single 
crystals with x = 1. Figure. Eta) shows a 0-scan through 
the (1 1 0) Bragg reflection measured with neutron 
diffraction using the SPINS spectrometer at the NIST 
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FIG. 2: (color online) (a) The bulk susceptibility, M/H, of an 
x — 1 single crystal as a function of temperature. Also plot- 
ted are results from powder samples. Inset: Low temperature 
specific heat of an x = 1 single crystal compared with powder 
data (with H = 0). (b) Plots of MT/H versus temperature, 
showing the anisotropy Xc > Xab at high temperatures. In- 
set: M versus H curves at 80 K and 300 K along different 
directions. 

Center for Neutron Research. The width of the scan of 
0.5°, which is resolution-limited, attests to the crystal 
homogeneity A high-resolution x-ray diffraction experi- 
ment on a small a; = 1.0 single crystal was performed on 
the ID20 beamline of the European Synchrotron Radia- 
tion Facility (ESRF). The scattering was performed in 
reflection geometry with the (10 1) reflection roughly 
perpendicular to the mount, with x-rays of energy of 
8.979 kcV (A = 1.381 A). The sample was cooled with 
a closed-cycle displex and was mounted on a four-circle 
goniometer. The measurements indicate that the sample 
remains in space group R3m down to the lowest measured 
temperature T ~ 2 K. The high x-ray flux available at a 
synchrotron is ideal to look for very subtle superlattice 
reflections that would arise if a VBS ground state resulted 
in a structural distortion. The 36-site VBS would lead 
to an enlarged unit cell that is a factor of VT2 longer on 
each side and rotated 90° from the original unit cell. We 
searched for the superlattice reflections along the high- 
symmetry (110) direction between the (4 4) and (5 1 4) 
peaks. This scan, shown in Fig.[3Jb), showed no observ- 
able superlattice peaks above the background. The only 
scattering features in this range were very weak powder 
peaks (roughly 250,000 times weaker than the strongest 
lattice reflections). These peaks were confirmed to arise 
from powder through f5-scans; they also displayed no tem- 
perature dependence and were somewhat broader than 
resolution. Contamination from these powder peaks, al- 
though exceptionally weak, is likely the limiting factor 
in setting an upper bound on the possibility of any su- 
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FIG. 3: (color online) (a) Neutron diffraction #-scan through 
the (110) Bragg reflection of an x = 1 single crystal measured 
with the SPINS neutron spectrometer, (b) Synchrotron x-ray 
diffraction intensity of a scan along the (4+e e 4) direction 
at three temperatures, 2 K, 20 K and 50 K. The extremely 
weak peaks observed in this range arise from powder contam- 
ination and have no temperature dependence. The lines are 
a simulation of the superlattice peaks arising from a 36-site 
valence bond solid with a 1% reduction in the bond length of 
the static dimer pairs, as discussed in the text. 

pcrlatticc reflections. The blue lines in Fig. [3Jb) are a 
calculation of the superlattice peaks expected if the cop- 
per ions making up the 6 static dimer pairs per super- 
cell (dimers that are in neither a pinwheel nor a perfect 
hexagon) were to move toward one another such that 
the Cu-Cu distance was reduced by 1%. It was assumed 
that the centering of the enlarged unit cells on adjacent 



kagomc planes is random, so that the superlattice peaks 
will actually be rods of scattering parallel to c. The lack 
of observed superlattice peaks at the expected positions 
indicates that any structural distortions due to this pro- 
posed supercell are below the 1% level. Note that the 
intensity of the superlattice peaks is proportional to the 
displacement squared in the limit of small displacement. 
Additional mesh scans (not shown) found no evidence of 
other superlattice peaks, as might arise from VBS states 
with different enlarged unit cells^. Of course, this ex- 
periment cannot rule out the presence of a VBS ground 
state that results in very little or no structural change. 

In summary, high quality single crystals of 
Zn x Cu4-a;(OH)6Cl2 have been synthesized and charac- 
terized. The bulk properties of x = 1 single crystals 
are consistent with the previously published powder 
results. The susceptibility measured along different 
crystallographic directions shows clear anisotropy. This 
indicates the presence of additional terms in the spin 
Hainiltonia.n, such small Dzyaloshinskii-Moriya 

interaction or exchange anisotropy. Synchrotron x-ray 
scattering experiments did not show evidence for the 
emergence of superlattice peaks at low temperatures. 
Hence, any possible lattice distortions associated with 
a valence bond solid are subtle, if they exist. Clearly, 
further measurements on these single crystal samples 
should help reveal the physics of the S = 1/2 kagome 
ground-state. 
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